* Re-Assessing Governmental Corruption of Religion: Where Political Theory Meets Empirical Evidence *
* Replication DO File *

*********************************
* All Christian Majority States *
*********************************
asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) reset

asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

**********************************************
* All Christian Majority States, Odds Ratios *
**********************************************
asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) reset

asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

*****************************
* Marginal Effects Analyses *
*****************************

meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

margins, dydx(lz06funding_l1) post

est store lz06_maj

meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

margins, dydx(lz06funding_l1) post

est store lz06_min

coefplot (lz06_maj, label(Majority) offset(-0.001)) (lz06_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Church) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(General Religious Funding) saving(lz06)



meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

margins, dydx(vfundingx_l1) post

est store vfund_maj

meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

margins, dydx(vfundingx_l1) post

est store vfund_min

coefplot (vfund_maj, label(Majority) offset(-0.001)) (vfund_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Church) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(Exclusive Religious Funding) saving(vfund)



meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

margins, dydx(lz07entanglement_l1) post

est store lz07_maj

meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

margins, dydx(lz07entanglement_l1) post

est store lz07_min

coefplot (lz07_maj, label(Majority) offset(-0.001)) (lz07_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Church) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(Religious Entanglement) saving(lz07)

gr combine lz06.gph vfund.gph lz07.gph, ycommon row(1) ysize(2) xsize(6) scale(1)
